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ABSTRACT 


This  paper  contains  an  analysis  of  the  errors  introduced  into  a  monopulse  radar 
system  by  employment  of  a  digital  matched  filter  containing  two  high  speed  analog-to- 
digital  converters  and  describes  errors  predicted  by  a  computer  simulation  of  the  sys¬ 
tem.  The  analysis  and  simulation  indicates  that  a  digital  monopulse  would  compare  very 
favorably  in  performance  with  an  analog  monopulse  system. 


This  technical  documentary  report  is  approved  for  distribution. 


Franklin  C.  Hudson,  Deputy  Chief 
Air  Force  Lincoln  Laboratory  Office 


The  digital  monopulse  radar  was  proposed  and  designed  by  E.  Gehrels,  G,  M.  Hyde, 
L.G.  Kraft,  and  J.  Arthur,  Group  31,  Lincoln  Laboratory. 
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1. 


INTRODUCTION 


A  block  diagram  of  the  Millstone  monopulse  radar  control  system  is  shown  in 
Fig.  1,  Azimuth  and  elevation  sum  and  difference  signals  (denoted  by  A  ,  A  , 

2^  ,  ,  and  )  of  frequency  1295mc  are  formed  at  the  antenna  from  the  four  re- 

ceived  signals.  The  four  sum  and  difference  frequencies  are  beat  down  to  28  me,  29mc, 
30mc,  and  31mc  respectively  near  the  receiving  antenna  and  then  linearly  combined  and 
transmitted  to  the  ground  through  a  single  channel.  The  doppler  is  subtracted  out  and 
the  signal  is  again  beat  down  by  mixing  it  with  a  30mc  signal.  The  sine  and  cosine  com¬ 
ponents  of  the  resulting  signal  are  then  separated  and  assuming  that  they  have  an  additive 
noise  component  the  resulting  two  signals  have  the  form  given  in  Eqs.  (1)  -  (2). 

y^  =  S  cos  27rf^t-Tsin  27rf^t+Ucos  27rf^t+V  sin27rf^t+Xcos  27rf2t+Ysin27rf2t+P+Noise  (1) 
y^  =  S  sin27rf^t+Tcos  27rf^t-Usin  27rf^t+V  cos  27rf^t-Xsin27rf2t+Y cos  27rf2t-Q+Noise  (2) 

where 

f^  =  10^  cycles/sec. 
f^  =  2  •  10^  cycles/sec. 

and  the  coordinate  pairs  (S,  T),  (U,  V),  (X,  Y),  and  (P,  Q)  are  the  complex  coordinates  of 
vectors  defined  as  follows. 

-  i  wt 

A  =Re[(S+iT)e  ]  =  S  sin  wt  +  Tcos  wt 
“  i  wt 

A  =Re[(U+iV)e  ]  =  Usinwt  +  Vcoswt 
A  Lt 

^Ortho  =  Xsinwt  +  Ycos  wt 

=Re[(P+iQ)e  =  Psin  wt  +  Q  cos  wt 
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A  digital  filter  is  employed  to  determine  S,  T,  U,  V,  X,  Y,  P,  and  Q  by  separat¬ 
ing  out  the  various  components  of  y  and  y  in  the  following  manner.  A  number,  m,  of 

^  1 

uniformally  spaced  samples  of  y^^  and  y^  one -fourth  of  a  usee  apart  (  —  of  a  cycle  with 
respect  to  the  Imc  components  of  y^  and  y^)  are  added  together  as  specified  by  Eqs. 

(3)  -  (10).  The  time  of  duration  of  the  sampled  pulse  is  nominally  2  millisec.  so  that  a 
total  of  8,000  samples  (m  =  8,000)  are  added  together  in  accordance  with  Eqs.  (3)  -  (10). 
The  resulting  sums  are  proportional  to  the  appropriate  components  of  y^^  and  y^  under 

^1  '^1  '^2 

4  •  2000 


1999 


n=0 


1999 


T  = 


n=0 


^2  )  '^1  ^)^)  -Y2  (^(2n+l)7rj  +y^  (^(2n+  |)7r 

4  •  2000 


(4) 


1999 


U  = 


n=0 


Yl  i(2n7r)  -y^ 


((2n+|)7r)  -y^ 


1  (2n+l)tr  I  Ty^  I  (2n+  — )7r  I  (5) 


4  •  2000 


1999 


V  - 


n=0 


^2  ■'■^1  '^2  (^^2n+l)7rj  -y^  ^(2n+|)7r^ 

4  •  2000 


(6) 


1999 


X  = 


n=0 


y^  (^2n7rj  -y^  ((2n+  ^)7rj  +y^  (^(2n+l)7rj  -  y^  ^(2n+  |)7r 

2  •  4000 


(7) 


1999 


Y  = 


y^  I  2n7r  j  -y 2  U2n+  2)'’^ j  +y 2  \  (2n+l)7ry 


■^2 


n=0 


2  •  4000 


(8) 
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1999 


P  =  y  Yj  \2nn  j+y^  \^(2n+  -)7r j  +y^  ^(2n+l)7r j  +y^  (^(2n+  -)7r^ 


n=0 


2  •  4000 


(9) 


1999 


Q  = 


n=0 


•y2  y2n7rj  -y^  (^(2n+ |)7rj  -  ((2n+l)7rj  -y^  (^(2n+  |)7rj  (10) 

2  •  4000 


(the  symbol  ~  indicates  the  estimate  of  the  signal) 

the  assumption  that  the  mean  of  m  samples  of  noise  is  zero.  For  example,  it  is  seen 
that  evaluation  of  the  expression  specified  in  Eq.  (3)  results  in  a  signal  proportional  to 
S  provided  that  the  mean  of  the  noise  samples  is  zero. 


S  = 


^  Yj  (^nTr)  +y^  ^(2n+ -y^  ^(2n+l)7r^ 


n=0 


Y2  ^(2n+2)\' 


4  •  2000 


(3) 


where 


2n.J 

=  S  +U  +  X 

+  P 

+  Noise 

(2n+  ^)7r^ 

=  S  -  U  -  Y 

-Q 

+  Noise 

(2n+l)7r  J 

=  S  +  U  -  X 

-  P 

-I-  Noise 

(2n+  |)7r  j 

=  S  -  U  + Y 

+  Q 

+  Noise 

7999 


(Noise). 


S  =  S  + 


i=0 


4  2000 
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The  antenna  azimuth  and  elevation  control  signals  may  be  obtained  from  the 
estimates  for  S,  T,  U,  V,  X,  and  Y  by  computing  the  normalized  projections  of  the 
vectors  specified  by  the  coordinate  pairs  (S,  T),  and  (U,  V)  on  the  vector  specified  by 
the  coordinate  pair  (X,  Y)  as  presented  in  Eqs.  (11)  -  (12).  A  subsequent  paper  will 
derive  these  equations  and  describe  the  computer  necessary  to  perform  these  calculations. 


Azimuth  signal  to  servo 


^AZ  ^  UX  +  VY 

y  9  9 

Ortho  X  +  Y 


(11) 


Elevation  signal  to  servo 


^EL  ^  SX  +  TY 

y  ~  9  9 

Ortho  X  +Y 


(12) 


2.  ERROR  ANALYSIS 

An  analysis  of  the  errors  in  the  estimates  of  the  amplitude  of  the  components, 

S,  T,  U,  V,  X,  Y,  P,  and  Q,  of  the  two  noisy  signals,  y^^  and  y^  is  now  examined. 

The  signals  are  of  the  form  given  in  Eqs.  (1)  -  (2).  Since  the  signals  are  sampled  syn- 

nTT 

chronously  at  a  time  when  27rf^t  =  (n  integral)  the  signals  can  be  considered  to  be 
of  the  form 

y(t.)  =  S  +  n(t.)  (13) 

th 

where  n(t.)  is  the  magnitude  of  the  noise  at  the  time  of  the  i  sample  and  S  is  the  com¬ 
ponent  of  y^  and  y^  that  is  to  be  estimated.  If  the  value  of  the  i*"^  quantized  sample  is 
denoted  by  x(t.),  then  the  estimate,  S,  of  the  noisy  signal  is 


4 


m 


i=l 


(14) 


Tlie  noise  is  assumed  to  be  additive  uncorrelated  from  sample  to  sample,  stationary  and 

2 

normally  distributed  with  zero  mean  and  variance,  a  . 

n 


n(t  =  (27ra^)‘^/^ 
i  /  n 


exp 


i  2,  ,  .  2 

-  -  n  (t  )  /  a 
Z  I  T\ 


(15) 


In  general,  a  quantizer  has  the  input-output  characteristic  shown  in  Fig.  2.  Thus, 
if  the  input  to  the  quantizer,  y(t.)  is  such  that  b.  ^  y(t^)  <  b^^^^  then  the  output  of  the  quan¬ 
tizer,  x(t  ),  is  a  where  b.  a.  ^  b.  .  ,  but  a.  is  not  necessarily  equal  to  ^i  ^iH-1  . 

11  1  1  iH-1  1  - - - 

It  should  be  noted  that  the  assumed  quantizer  is  symmetric: 

a  .  =  -a 
-1  1 

and 

b  .  =  -b., 

“1  1 


but  (a.  -a.)  is  not  necessarily  equal  to  (a. -a.)  nor  is  (b.  -b.)  necessarily  equal  to 

i+k  1  j-Hc  j  i+k  1 

(b.  ,  “b.)  (i.e.,  equal  quantizing  levels  are  not  necessarily  assumed).  The  quantizer 

j 

input-output  characteristic  is  completely  determined  by  the  two  column  vectors 


a. 

0 

1 

a^ 

b^ 

2 

2 

and  B  = 

Tlie  vector  A  determines  the 't  output  levels  and  in  order  to  be  easily  realizable  it  must 
be  of  the  form 


5 


A  =  a 


1 

2 

3 


I 


(equal  quantizing  steps)  or  of  the  form 


where  L  is  an  integer  (logarithmic  quantizing  steps). 

The  vector  B  determines  the  division  of  the  input  into  various  levels,  and  the  components 
of  B  may  take  on  any  values  without  presenting  significant  implementation  problems. 
(These  levels  are  usually  determined  by  the  value  of  resistors). 

The  probability  density  of  the  quantizer  output,  x,  is 


where 


k=— t 


=  Prob.  = 


Prob.  Ib^  <  y  <  b^^jJ 


b 


k+1 


-1/2 


b, 


k 


n 


It  is  convenient  to  define  the  following  column  vectors: 
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— 

Pi 

P-i 

1 

^1 

2 

P2 

Q  = 

P-2 

and  = 

^2 

• 

■  2 

Pf. 

_P-t 

The  average  value  of  the  quantizer  output  E[x]  is  then 


I 

=  I  Vk 
1^=-^ 

=  aV-Q) 


(16) 


and  the  variance  ct  is 

X 


-t 


k=-T 


=  aJ(p-k5)  -  [aV-Q)]^ 


(17) 


where  A  is  the  transpose  of  A,  the  vector  defining  the  quantizer  output  levels  and  P 
and  Q  are  functions  of  the  vector  B.  Since  the  output  of  the  sampler  is  independent 
(uncorrelated)  from  sample  to  sample,  the  output  of  the  quantizer  is  also  independent 
from  sample  to  sample. 

The  sum, 

S  = 

i=l 


m 


1 

m 


x(t.) 


(14) 


stored  in  the  accumulator  is  used  as  an  estimate  of  S.  The  average  value  of  S  is 


m 


E[S]  =  E 


—  >  x(t.) 

m  /_j  i 

i=l 


m 

=  ^y  E[x]  =  E[x] 
m 

i=l 
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Therefore  the  average  value  of  the  estimate  of  S  is 


E{S]  =  A  (P-Q). 


(18) 


The  variance  of  S  is 


0-2  =  E  [s^ 
S 


■  1  m  m 

'  T 

E 

— 2  .S  Z)  x(t.)x(t.) 

_  m  1=1  J=1  1  j 

A  (P-Q) 

However 


E[x(t.)]  E[x(0],  i5^j 


E[x(t.)x(0]  =  ^ 


E(x  (t.)] 


1  =  J 


(  T  2 

[A  (P-Q)]  ,  i-^j 


A2(P-+Q),  i  =  j- 


Therefore 

2 

0-^ 

S 


1 

— 2 
m 


mAj(P+Q)  +  (m^-m)  [a'^(P-Q)]^ 


T  2 

-  (A  (P-Q)] 


J_ 

m 


aJ  (P-+Q)  -  [aV-Q)]^ 


(19) 


1  2 

—  a  . 
m  X 


The  relative  error,  r,’  is  defined  as 
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r  = 


(20) 


s  -  s 
s 


The  average  value  of  this  error,  E  [  r] ,  is 


E[r] 


S  -  E[?] 

S 


^  S  -  A^(P-Q) 

S 

and  the  variance  of  the  relative  error  is 

2  2  2 
=  E[r  ]  -  E  [r] 


CT  2 

S' 


It  is  seen  from  Eq.  (21)  that  the  average  value  of  the  relative  error  resulting  from  em¬ 
ployment  of  a  digital  filter  is  independent  of  the  number  of  samples  and  dependent  only 
on  the  characteristics  of  the  quantizer.  On  the  other  hand,  the  variance  of  the  relative 
error  (Eq.  (22)  )  is  inversely  proportional  to  the  number  of  samples,  m,  and  depends  on 
both  the  number  of  samples  and  the  characteristics  of  the  quantizer.  Therefore,  increas¬ 
ing  m,  the  number  of  samples,  will  decrease  the  variance  of  the  relative  error,  but  will 
have  no  effect  upon  the  average  value  of  the  relative  error. 
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For  large  m,  the  Central  Limit  Theorem  may  be  invoked.  Then  S  is  normally 


distributed  with  variance  and  mean  E[S]  . 

S 


Thus  for  large  m 


p(S)  =  I  27ra 


■1/2 


exp  < 


1^ 

2 


(S  -  E[S]) 


a 


(23) 


and 


P(r)  = 


r  /  “1  Y 

-  1/2 

1  (r-E(r))^  1 

exp  ^ 

"2  2  f 

(24) 


The  square  root  of  the  average  value  of  the  squared  relative  error^E(r  ),  can 
be  used  as  a  suitable  measure  of  the  error  in  estimating  S.  A  computer  program  to 
evaluate  E(S),  E(r),  a  ,  andl^ (r^)  =VE^(r)  -f  cr^  as  specified  in  Eqs.  (18)  -  (22)  for  any 

I*  ^  jjc 

set  of  values  of  m,  a  ,  S,  A,  and  B  has  been  prepared. 

n 

Consideration  of  data  obtained  through  employment  of  the  computer  program 
resulted  in  the  following  observations. 

1.  Consideration  should  be  given  to  utilization  of  an  A/D  converter  having  loga¬ 
rithmically  spaced  outputs  in  radar  systems  of  this  type.  For  example,  it  was  found 

that  for  ranges  of  m,  a  ,  S,  A,  and  B  pertinent  to  this  radar  system  smaller  values  of 
2  ^ 

E(r  )  result  if  the  A/D  converter  employs  7  logarithmically  spaced  quantization  levels 
(7  bits)  rather  than  32  equally  spaced  output  levels  (5  bits)  because  of  the  greater  dynamic 
range  of  the  logarithmically  spaced  output  A/D  converter.  It  should  be  noted  that  at  the 
4  me  sampling  rate  employed  a  7  level  logarithmically  spaced  output  A/D  converter  is 
probably  simpler  to  construct  than  the  latter  5  bit  32  level  device. 

2.  Reducing  the  number  of  bits  employed  in  the  A/D  converter  to  less  than  5 
bits  in  order  to  make  feasible  the  utilization  of  a  higher  sampling  frequency  (at  present 

* 

Acknowledgement  is  made  to  R.  Carroll  for  preparation  of  these  programs 
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It  is  very  difficult  to  construct  A/D  converters  of  greater  than  4  me  sampling  frequency 

which  contain  large  numbers  of  quantization  levels)  results  in  an  increase  in  the  value 
2 

of  E(r  )  for  the  range  of  m,  S,  A,  and  B  relevant  to  this  radar  system.  The  effect 
of  decreasing  the  number  of  quantization  levels  and  simultaneously  increasing  the  sam¬ 
pling  frequency  (and  thus  increasing  m  by  beating  down  the  1295mc  signals  to  frequencies 

2 

greater  than  Imc  and  2mc)  is  to  increase  the  component  of  E(r  )  due  to  the  nonzero 

2 

average  of  the  relative  error  E(r)  and  decrease  the  component  of  E(r  )  arising  from  a  . 

2  ^ 

3.  As  expected,  the  component  of  E(r  )  arising  from  the  nonzero  value  of  E(r) 

can  be  made  relatively  small  near  values  of  S  corresponding  to  quantization  levels 

(components  of  A) 

It  should  be  noted  that  this  error  analysis  presents  the  errors  in  S,  T,  U,  V,  X, 
Y,  P,  and  Q  rather  than  the  expected  errors  in  the  azimuth  and  elevation  signals  to  the 
servo  system.  In  order  to  obtain  this  latter  information,  and  in  order  to  take  into 
account  other  factors  contributing  to  errors  in  aiming  the  radar  antenna  such  as  not 
completely  compensating  for  Doppler  frequency  shifts  in  the  received  signals,  a  com¬ 
puter  simulation  of  the  system  was  programmed.  The  results  of  the  simulation  indicated 
that  if  Raytheon  4mc  7  bit  A/D  converters  were  employed  then  the  radar  system  inaccu¬ 
racies  resulting  from  Doppler  errors  (resulting  in  an  incorrect  sampling  frequency) 
would  be  a  significant  portion  of  the  overall  system  inaccuracies  and  the  expected  errors 
in  computing  the  servo  signals  would  be  a  few  percent  for  the  Doppler  errors  and  signal - 
to-noise  ratios  present  in  the  radar  system. 
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